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Abstract. We carry out a general study of the stability of astrophysical flows that appear steady in a uniformly 
rotating frame. Such a flow might correspond to a stellar pulsation mode or an accretion disk with a free global 
distortion giving it finite eccentricity . 

We consider perturbations arbitrarily localized in the neighbourhood of unperturbed fluid streamlines. When 
conditions do not vary around them, perturbations take the form of oscillatory inertial or gravity modes. However, 
when conditions do vary so that a circulating fluid element is subject to periodic variations, parametric instability 
may occur. For nearly circular streamlines, the dense spectra associated with inertial or gravity modes ensure that 
resonance conditions can always be satisfled when twice the period of circulation round a streamline falls within. 
We apply our formalism to a differentially rotating disk for which the streamlines are Keplerian ellipses, with free 
eccentricity up to 0.7, which do not precess in an inertial frame. We show that for small e, the instability involves 
parametric excitation of two modes with azimuthal mode number differing by unity in magnitude which have a 
period of twice the period of variation as viewed from a circulating unperturbed fluid element. Instability persists 
over a widening range of wave numbers with increasing growth rates for larger eccentricities. 
The nonlinear outcome is studied in a follow up paper which indicates development of small scale subsonic 
turbulence. 
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1. Introduction 

In Astrophysics one is often concerned with the stability of 
a system undergoing some regular global steady motion. 
Possible examples are stars undergoing a pulsation that 
may be either freely excited or forced through tidal inter- 
action in a binary system. Other examples, which are more 
the focus of this paper, are differentially rotating disks 
which have been either set up or perturbed into a state in 
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which the orbital motion about the centre is non circular. 
Systems of this type are known to be susceptible to para- 
metric instabilities which involve a resonance between the 
periods of normal modes of oscillation and a period associ- 
ated with the unperturbed motion. For examples of stud- 
ies related to stellar oscillations see Vandakurov(1967), 
Papaloizou (1973), Papaloizou & Pringle 1980), Kumar & 
Goodman (1996), and Wu & Goldreich (2001). 
The presence of accretion disks in which the orbits are ec- 
centric have been suggested in a number of contexts. The 
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production of such disks by the tidal disruption of stars or- 
biting around AGN has been considered by Gurzadyan & 
Ozernoy (1979). A model to explain the long term spectral 
variability of Be stars as being due to global one-armed 
(to = 1) oscillations of equatorial disks has been devel- 
oped (eg. Okazaki 1991; Papaloizou, Savonije & Henrichs 

1992) . Eccentric disks have also been proposed as an ex- 
planation for the superhump phenomenon in SU UMA 
stars (Whitehurst 1988). 

A general theory to describe the slow viscous evolution of 
disks with modest to large eccentricities has been devel- 
oped by Ogilvie (2001). This confirms the possibility of 
long lived structures of this type within that context. 
The stability of the non circular motions induced in an 
accretion disk by a binary companion has been studied in 
pioneering work by Goodman (1993) and Ryu & Goodman 
(1994). These authors studied the fast to = 2 forced dis- 
tortions that propagate around with the companion. A 
parametric instability was found that lead to turbulence 
in the context of local two dimensional modelling in a 
shearing box. 

In this paper we are concerned with extending studies of 
the stability of disks with non circular motions to more 
general situations such as global disks with free global 
TO = 1 modes which are almost stationary in an inertial 
frame. 

The recent discovery of a number of extrasolar giant plan- 
ets orbiting around nearby solar-type stars ( Mayor & 
Queloz 1995; Marcy & Butler 1998, 2000) with high orbital 
eccentricities in the range 0.0 < e < 0.93 has led to the 
consideration of the proposal that a giant planet orbiting 
in and interacting with a protostellar disk (eg. Goldreich & 
Tremaine 1980; Ward 1997) could undergo an instability 
that would cause its eccentricity to increase. (Artymowicz 
1992; Lin & Papaloizou 1993; Papaloizou Nelson & Masset 
2001; Goldreich & Sari 2003). However, the free modes in 
such an interaction are likely to be coupled to the disk 
which participates with significant eccentricity in a global 
TO = 1 mode. (Papaloizou 2002, Goldreich & Sari 2003). 
Accordingly the behaviour and in particular the damp- 
ing of free eccentricity in differentially rotating disks be- 
comes an issue. It is one of the purposes of this and a 
follow up paper (Papaloizou 2004) to study the instabil- 
ities of free global disk normal modes corresponding to 
free eccentricity and their non linear outcome. Although 
there may be major regions of protostellar disks that sup- 
port magnetic fields leading to turbulence through the 
Magnetorotational instability (Balbus & Hawley 1991), 
because of dependence on the existence of external sources 
of ionization and complex chemistry, the extent of such re- 
gions is highly uncertain ( eg. Gammie 1996 ; Promang, 
Terquem & Balbus 2002). Accordingly for this first study 
we have considered disks with eccentricity in the hydrody- 
namic limit. This is of general theoretical interest because 
even in that limit there are potential local instabilities of 
disks with non circular but regular streamlines that could 
lead to turbulence and enhanced dissipation (Goodman 

1993) . Although uncertain, it is likely that the presence 



of magnetic fields would enhance any decay rate found in 
the hydrodynamic limit. However, this will be studied in 
a separate investigation. 

In this paper we develop a general formalism for the study 
of the local stability of astrophysical flows possibly in stars 
or disks that appear steady in a uniformly rotating frame. 
We adopt a local Lagrangian approach and consider per- 
turbations localized near unperturbed fluid streamlines. 
These are shown to be governed by a pair of second order 
ordinary differential equations that give the evolution of 
the localized disturbance in the neighbourhood of a fluid 
element as a function of time. For streamlines on which 
conditions vary as a fluid element circulates, the coeffi- 
cients are periodic functions of time. The parametric ex- 
citation of inertial or gravity modes readily occurs. 
We solve these equations and delineate the instability for 
an example for which the streamlines are Keplerian el- 
lipses, with free eccentricity up to 0.7. The nonlinear out- 
come of such instabilities in disks supporting global modes 
which endow them with eccentricity is studied in a follow 
up paper which indicates development of small scale sub- 
sonic turbulence. 

The plan of the paper is as follows: In section 2 we give 
the Basic equations. In section 3 we describe the steady 
state flow conditions of the systems whose stability we 
analyse, in the context of a model with streamlines that 
are Keplerian ellipses. In section 4 we discuss the general 
linear stability analysis and go on to describe the solution 
procedure used in the local limit in general terms relat- 
ing the spectra found in the local problem to that of the 
global one in appendix 3. In section 5 this is applied to 
the case when there are circular streamlines in the un- 
perturbed steady state. In section 6 this is extended to 
a steady state model with streamlines taking the form 
of Keplerian ellipses. Numerical and analytic solution of 
the governing equations are described in section 7 for the 
Keplerian streamline case while emphasizing the general 
nature of the behaviour. Finally in section 8 we discuss 
our results 

2. Basic Equations 

We work in a rotating frame which rotates with angular 

velocity ftp relative to an inertial frame. 

We adopt a system of cylindrical coordinates {r,ip,z) 

for which the origin is at the centre of mass and the 
angular velocity is directed along the z axis such that 
rjp = (0,0, f7p). 

In the absence of magnetic fields, the equation of motion 
may be written in the form 

dv VP 
—+v-Vv+2ftpXw+ftpx{np X r) = F = V'J'.(l) 

Here the velocity is v = (u^, f,^, i',~). The left hand side 
includes contributions from the acceleration of a fluid ele- 
ment and coriolis and centrifugal accelerations arising be- 
cause of the rotating frame. The acceleration of a fluid 
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element may be written using the convcctive derivative 
taken moving with the fluid in the form 



Dv 

'm 



dv 
'di 



+ V • Vv. 



(2) 



The system is subject to a force per unit mass F which 
has components due to pressure. P. and gravity, with the 
gravitational potential being ^. The density is p. 
In addition we have the continuity equation 

dp 



dt 



+ V • (pv) = 0. 



(3) 



3. Steady State Configurations 



We here consider the stability of steady state configura- 
tions for which d/dt = 0. For these equation (1) gives 



V • Vv + 2flp xv + r2px(r2pxr) 



VP 

~ V*.(4) 



In general this allows for a steady state flow v as viewed 
in a rotating frame. 

A very simple possible motion is non uniform rotation for 
which V = {0,v^{r),0) corresponding to circular motion 
about the z axis at a rate which is a function of r. 
Other examples could correspond to a normal mode of 
stellar oscillation as viewed in a frame rotating with the 
modal pattern speed or an accretion disk with a global 
eccentricity. 

3.1. Cylindrical Disks 

In the context of differentially rotating accretion disks we 
consider fluid circulating in Keplerian ellipses in a cylin- 
drical potential. For this Hp = and \E' = — GM*/r. Here 
M* is a central mass and G the gravitational constant. 
Disk self gravity is neglected. If in addition the pressure, 
P, is constant, we have 



Vv = -V*. 



(5) 



This model can approximate conditions near the midplane 
of a vertically stratified disk if dvz/dz = there and with 
a slight modification if it is not. 

Equation (5) applies to particle motion under only the po- 
tential Thus a solution can be found from the Keplerian 
trajectory 



a(l-e2) 



1 + e cos ip, 



(6) 



where the semi major axis a and the eccentricity e are 
constant on a streamline being a particle trajectory. On 
such a trajectory the velocity components are given by 



Vr = esmifA 



and 



r(l -|- ecosf^) 



IGM^{1 + ecosip) 



(7) 



(8) 



The eccentricity is constant on a stream line and may be 
specified as an arbitrary function of a. Equation (6) may 
then be solved for e as a function of r and (p. The density 
may then be found from the steady state continuity equa- 
tion V • (pv) = 0. In the special case when e is constant, 
the continuity equation is satisfied for p being a general 
function of its arguments of the form 



p = p (r(l -|- e cos ip),z) . 



(9) 



For the work presented later on we adopt a locally isother- 
mal equation of state such that 



P = pclir). 



(10) 



The sound speed Cg (r) is taken to be a flxed function of r. 
As above we adopt a cylindrical disk model for which verti- 
cal stratification is neglected (eg. Hawley 2000; Steinacker 
& Papaloizou 2002) and the potential is independent of z 
such that $ = —GM^,/r. We note again that this model 
can describe conditions close to the midplane of a stratified 
disk which is all that is required in order to perform a local 
stability analysis of the type we consider here. When p is 
constant and cl{r) oc 1/r, Keplerian ellipses are possible 
streamlines that form a stationary pattern in the inertial 
frame. 

3.2. Properties of Streamlines 

We study the stability of the steady state configurations 
described above, below. We focus on perturbations that 
are localized on particular streamlines and we assume that 
the configuration of steady state streamlines is smooth 
enough that at least locally the fluid volume can be filled 
with stream tubes that can be used to define a volume 
element dV = dAds, with dA being the stream tube area 
element normal to a streamline with element of length ds. 
Then the steady flow condition V • pv = on a stream- 
line may be written d{p\v\A)/ds = 0. We also suppose 
that two independent invariants can be defined on each 
streamline. For such invariants Q,i = 1,2 



V • VCi = 0. 



(11) 



In the simple case of differential rotation about the z axis, 
the stream lines are circles centred on the z axis. The 
invariants can be taken to be r and z, the radius of a 
particular circle and its height z respectively. Similarly 
for the cylindrical potential models with streamlines being 
Keplerian ellipses. The invariants could be taken to be z 
along with the semi-major axis or eccentricity, both of 
which are constant on a streamline. 

4. Linear Stability Analysis 

4.1. General Formalism 

We now study the stability of steady state flows of the 
type introduced above. The most convenient formalism 
is the Lagrangian approach developed by Lynden-Bell & 
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Ostriker (1967). Following these authors we introduce the 
Lagrangian variation A such that any state variable Q is 
perturbed following a fluid element from its value appro- 
priate to the steady state such that 



Q^Q + AQ. 



(12) 



The Lagrangian displacement is introduced as the 
Lagrangian variation of the position vector of a fluid ele- 
ment so that 



Ar = ^ 

The Lagrangian variation of the velocity is then 

where the convective derivative of ^ is 
Dt dt ^ 

and as shown by Lynden-Bell and Ostriker(1967) 
£>Av 



Dt ' 



(13) 



(14) 



(15) 



(16) 



The Eulerian variation, Q' in Q seen in a fixed coordinate 
system resulting from AQ is given by 



Q' = AQ - ^ • VQ. 



(17) 



To obtain the basic governing equations for the linear sta- 
bility problem, we take the Lagrangian variation of the 
equation of motion (1) which gives. 



_| + X + fip X (Jlp X ^) = AF. 



Using (17) we obtain 

VP' p' 



(18) 



AF 



/VP 

^^VP-^- V f — + v* 



We comment that because of the local nature of the per- 
turbations considered, we use the Cowling approximation 
in which variation of the gravitational potential is ne- 
glected. 

The Lagrangian variation in the density is related to the 
Lagrangian displacement through 



C • Vp = -pV • 



(20) 



The pressure and density variations could be related by 
the adiabatic condition 



AP _ TAp 



(21) 



where F is the first adiabatic exponent. For simplicity we 
shall assume this to be constant below, however this is not 
an essential requirement for what follows. 
One can seek solutions for which the time dependence is 
through a factor exp{iat) where a is the eigenfrequency 



associated with a normal mode and after taking all time 
derivatives divide this factor out. Then 



and 



-a^ + 2iav • V| V • V (v • V|)) . 



(22) 



(23) 



Adopting equations (22 - 23), the linear stability problem 
(18) can be reduced to the operator equation (see Lynden- 
BeU & Ostriker 1967) 



L(|) - -^2^ + <tB(|) -f C(|) = 0, 



(24) 



where, if the adiabatic condition (21) is used, the operators 
B, C are self adjoint with weight p. That is for two arbi- 
trary displacements (^, 77) and standard boundary condi- 
tions 



(77,C(0)= / pv*-C{$)dV : 



Jv 



C*{r,*)dV, (25) 



with the integration being taken over the domain of the 
fluid and a similar relation holding for B. This also defines 
the inner product when C — I, the identity. 
This self adjoint property, although useful for some pur- 
poses is not essential to the analysis carried out in this 
paper which can in fact be extended to equations of state 
for which it does not hold (see below). Accordingly we 
shall not consider the self-adjoint formalism further here. 

4.2. Local Analysis 

We here consider perturbations that are localized on 
streamlines. To do this we suppose that any perturbation 
quantity is of the form 



(19) AQ = (AO)oH/(Ci,C2)cxp(iA5(Ci,C2)). 



(26) 



Here the phase function S and the positive definite local- 
ization function W are functions of the streamline invari- 
ants Ci and C,2- The constant A is a large parameter. The 
effective wave number is 



k = AVS'. 



(27) 



On a particular streamline with Ci = Cio and (2 = C20 
say, this is a general linear combination of vectors in the 
VCi and VC2 directions. These are independent because 
Ci and C2 are. The parameter A can be used to make the 
length scale l/|k| ex 1/A arbitrarily small. The localiza- 
tion function W is chosen to approach zero as one moves 
away from a particular streamline on a length scale that 
approaches zero as A ^ 00 but more slowly than l/|k| say 
as oc l/Ai/2. 

The amplitude factor (A(5)o is assumed to be a function of 
time and slowly varying in space. Because of the localiza- 
tion this is eff'ectively along the streamline with ^1 = Cio 
and C2 = C20- 
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Adopting the above form of perturbations, because of the 
rapid variation of the complex phase XS, one can make the 
standard assumption of local analysis that only the vari- 
ation of this needs to be considered when taking spatial 
derivatives. The only exception occurs when the variation 
is only along streamlines. This happens when the deriva- 
tives in an expression occur in the form v • V as in for 
example equations(22 - 23). Then, because v • \/S = 0, 
the contribution v • V(AQo) must be retained. 
Proceeding in this way one may obtain a set of govern- 
ing equations for quantities on a streamline on which the 
perturbation is localized. Because the exponential and lo- 
calization factors cancel out we drop the subscript zero 
from the perturbation amplitudes. Apart from such a jus- 
tification as in standard local analysis one can also prove 
that the eigenvalues a one obtains belong to the spectrum 
of the global linear stability problem (see appendix 3 and 
also Papaloizou & Pringle 1982; Terquem & Papaloizou 
1996 for a discussion of related problems). 

4.3. Reduction of the Linearized Equations 

Before embarking on a local analysis, we perform a prior 
reduction of the stability problem to a convenient form 
without approximation. 
We note that 



^+20px-^=a. 
a= AF-fip X (J7p X $,) 



(28) 
(29) 



Using equations (4) and (19 - 21), after some algebra, we 
can express a in the form 



a = — v(4V^rf^-^l 



p \hJ p ^ \rp p J 

-hi- V(v- Vv-t-2f2p X v), (30) 
where H = P^/'^ . 

Using equations (19 - 21) we may also write 



TP 



(31) 



In what follows we find it convenient to decompose 
into solenoidal and irrotational parts such that 



^ Vx A VcP 



(32) 



The powers of the large parameter A are introduced so as 
to anticipate the order of terms of different types. Thus 
it will turn out that is mainly solenoidal with an ir- 
rotational correction smaller in magnitude by a factor of 
order A~^. Prom equation(31) it then follows that 



P' = 
and 



rp\7^(, 

HX3 



(33) 



a=^V 



fTPV^(j,\ 




/VP 




\ ) 


+ — 1- 
p 


Vpp 





+C • V(v • Vv + 2J2p X v). 



4.3.1. Effect of the Equation of State 



(34) 



We comment here that the above analysis assumed an 
adiabatic condition for which the stability problem has a 
self-adjoint character (eg. Lynden-Bell & Ostriker 1967). 
However, because of the local nature of the formalism con- 
sidered here we may also consider other equations of state 
for which the formal self-adjoint character does not hold 
in the fully global problem while it is recovered in the local 
problem (see section 7). 

In particular, with a later application to disk configura- 
tions in mind we consider equations of state of the form 
P = pc^ (r) . Here c is the sound speed which is specified to 
be an arbitrary function of r. The equation of state also 
holds for perturbations so that P' = p'c^. 
Then in place of equation (31) we have 



P' = -c'V ■ ipC). 



(35) 



This is of the same form if we identify H = p, and 
TP = He?. With this identification we now get instead 

of equation (30) that 



■I- V(v- Vv-|-2np X v). (36) 



This differs by the positioning of and the the absence 
of the second entropy gradient term. However, the posi- 
tioning of is immaterial for the local analysis presented 
(it may be taken through the differentiation) and thus the 
formalism is the same as in the adiabatic equation of state 
case when the unperturbed entropy is constant or equiva- 
lently VP/(rP) = Vp/p. 

4.4. Solution Procedure in the Local Limit 

Although we introduced the small parameter A~^, above 
we have not yet used this smallness to simplify the sta- 
bility problem. We here reduce the equation (28) together 
with equation (34) to leading order in A~^, which is here 
equivalent to taking the limit A""*^ 0. 
We begin by noting that to leading order in A~^, is 
solenoidal ( V • {^H) = 0) as follows from considering the 
right hand side of (32). 

To specify a form for that is solenoidal to leading order 
in A^^, we note that to this order only variations of S need 
to be considered (see section 4.2) so that this is equivalent 
to requiring that k • | = 0. 

In this case | may be written as a linear combination of 
two linearly independent unit vectors bi and b2 that are 
perpendicular to k oc VS. 
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To specify two such unit vectors we first use the fact that 
V • \/S = and adopt 



bi = 



We can then adopt 

V X VS* 



|v X V^l 



(37) 



(38) 



which gives a unit vector orthogonal to both v and k. 
In adopting the above basis we make the assumption that 
there are no stagnation points in the flow at which |v| = 0. 
In this case may then adopt the decomposition 



^ = abi+/Jb2, 



(39) 



where a and /? are of the general form given by equation 
(26). 

The solution strategy is then to substitute the form of the 
displacement (39) into equation (28) together with equa- 
tion (34), to take the scalar product of the result with 
respect to bi and b2 respectively. This results in two equa- 
tions to determine a and /3 on the streamline on which the 
perturbation is localized. 

Noting that as the dependence on and ^2 can be eventu- 
ally factored out and dropped and only convective deriva- 
tives are involved, in a Lagrangian description we can ulti- 
mately regard a and /3 as functions only of time on a fixed 
streamline. Accordingly, these form a pair of coupled or- 
dinary differential equations. 

We begin by applying this procedure ro the left hand side 
of equation (28) which yields after straightforward alge- 
bra: 



a + 2CoP + Cna + Ci2/3 = a • bi 
and 



(40) 



/3-2CoQ! + C2ia + C22/3 = a-b2. (41) 

Here in order to condense the notation, we set for the first 
and second order convective derivatives operating on any 
quantity Q 

DQ 
Dt 

and 



= Q 



(42) 



= Q. 



(43) 



It is particularly useful to recall that for time independent 
quantities such as the unit vectors introduced above 



DO 



(44) 



We further have 

Co = b2-bi-h(ripxb2)-bi = -brb2-(nf xbi)-b2,(45) 

and for i = 1,2 and j = 1, 2 

dj = hj ■ hi + {2np X i),) • bi. (46) 



Using equation( 34) we can find a • bi, and a • b2 which 
occur on the right hand sides of equations (40) and (40) 
respectively. 

We begin by noting that in the limit A""'^ ^ 0, V reduces 
to multiplication by iX\/S. Thus in this limit 



a = — I 



/rp|V5pvs' 



VP, fVP Vp 



\ PH 
+i ■ V(v- Vv-|-2np X v), 

with the form of ^ being given by equation (39). 
Using the above we find 

a-bi = £>iia; + -D12/3 
and 

a • b2 = £'210; + £>22/3, 
where for i = 1,2 and i = 1,2 

Dij = hi ■ {hj ■ V(v • Vv + 2np X v)) 



(47) 



(48) 



(49) 



+hi 



VP, 



VP 
TP 



P 



(50) 



Thus the governing equations for a and /3 may be cast in 
the final form 



a + 2CQl3 = Eixa + E-^20 
and 

'l3-2CQa = E2ia + E22l3, 



(51) 



(52) 



where Ei,j = Di,j — Cij. 

To complete the solution the scalar (j) may be deter- 
mined by taking the scalar product of equation(28) with 
k. However, this is not needed for our purposes so we do 
not consider it further here. 

We comment that the following identities, valid for any 
unit vectors, are useful for evaluating the coefficients in 
the governing equations for a and /3 that involve second 
time derivatives: 



b, -b, = -|bi|^ 
and 

hi ■ hj = —hi ■ hj + 



Djhi ■ hj) 
Dt 



(53) 



(54) 



The governing equations (51- 52) consist of a pair of sec- 
ond order differential equations for the amplitudes {a, (3) 
as a function of time. As the above is a Lagrangian descrip- 
tion the coefficients appearing in these equations generi- 
cally are periodic functions of time. This is because these 
coefficients describe fluid effects as seen by a ffuid element 
as it circulates periodically on an unperturbed streamline. 
Only when these conditions are invariant such as on cir- 
cular streamlines in a system with cylindrical symmetry 
are the coefficients constant. In the generic case the pe- 
riodicity of the coefficients leads to a system similar to 
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the well known cases of Matliicu's or Hill's equation (eg. 
Whittaker & Watson 1996). Accordingly we expect to en- 
counter parametric instability when a normal mode fre- 
quency is commensurable with the orbital frequency as- 
sociated with circulating around a streamline. The low- 
est order relation is when the normal mode frequency is 
one half of the orbital frequency. This is easy to demon- 
strate when the system is almost cylindrically symmetric 
and the streamlines depart from circles by a small amount 
(see below and appendix 2). Furthermore commcnsurabil- 
ity conditions can always be satisfied when the underly- 
ing modes are gravity /inertial modes because, as is well 
known, their spectrum is dense (eg. Papaloizou & Pringlc 
1982; Lin Papaloizou & Kley 1993). This density comes 
about because in the local limit the oscillation frequencies 
depend only on wavcnumbcr ratios so that wavcnumbcr 
components may be allowed to tend to infinity at fixed 
ratios while not affecting the oscillation frequency. One 
can use that property to relate the spectrum of the lo- 
cal problem to that of the global problem ( see appendix 
3). We illustrate the above remarks by considering a sys- 
tem with cylindrical symmetry and a system for which the 
streamlines are Keplerian ellipses of arbitrary eccentricity. 
Of course in the latter case when the eccentricity is small 
the streamlines deviate slightly from circles. 

5. The Case of Circular Streamlines in the Steady 
State 

In this case, Clp = 0, and the invariants on a streamline 

arc taken to be (i = r. and (^2 = z- Because the orbits are 
circles we have r = a, where a is the semi-major axis. 
For the phase function we adopt S = KrT + K^z with K,. 

and Kz being constants. 

This gives the wavenumber k = X\/S = {XKr, 0, XK^). 
For this steady state model v = {0,v^{r),0), equation 
(37) gives hi = ip being the unit vector in the azimuthal 
direction. Then from equation (38) we obtain ba = {K^r— 
Kri.) I K, where f and k are unit vectors in the radial and 
vertical directions respectively and K = \J K"^ + K^. 
Having obtained these unit vectors the forms of equations 
(51) and (52) can be found by direct calculation. 
Using Dr/Dt = v • Vf = v,p(f/r, and D(p/Dt = v • \I(p = 



—v^r/r, we find setting Q. 
streamlines, that 



'W 
and 

Dt 



K ■ 



V(p/r which is constant on 



(55) 



(56) 



From these it is easy to obtain Co = bi • b2 = K^^/K, 
bi • b2 = 1)2 • bi = 0, bi • bi = — f2^, and b2 • b2 = 

Further we find that for | given by equation (39) that 



and 



b2 • (^ • V(v • Vv)) = - 



(58) 



In addition for i 1, 2, j = 1, 2, 



VP, /VP 



P " VTP 
• (k X VP) • 



Vp 
P 



|k|p 



(59) 



with 5ij denoting the Kronnecker delta. 
Using (57-59) and above, equations (51) and (52) are re- 
duced to the simple forms 



K 

P = - 7pr^ + ^ ] P, 



K 



where 



k X 



(60) 
(61) 



In this case because there is no variation of state variables 
around streamlines the governing equations have constant 
coefficients and we may look for normal mode solutions 
for which a and (3 are a exp{iat) where a is the normal 
mode frequency. 

Setting a = ao exp(icri), /3 = /3oexp(icrf) equations (60- 
62) give a non trivial solution such that 



a: = — — p 



Ka 



and 



^2 _ 2n{rn' 



2n)Kl 



(63) 



(64) 



bi • (C • V(v • Vv)) 



(57) 



This gives the well known local dispersion for iner- 

tial/gravity modes in a differentially rotating system (see 
eg. Tassoul 1978 and references therein for a discussion). 
There are two contributions. The first is from the rota- 
tion and depends on the gradient of angular momentum 
which must be positive for stability (Rayleigh's criterion) 
and the second is from the entropy gradient and leads to 
gravity waves when stable {v"^ > 0) or convection when 
convectively unstable (Schwarzchild's criterion ). We shall 
always assume that the system is stable to convection and 
satisfies Rayleigh's criterion. Accordingly we consider sys- 
tems that are dynamically stable in this sense. 
However, as we shall see below parametric instability can 
occiir when state variables vary around streamlines mak- 
ing a perturbed fluid element subject to periodic restoring 
forces even though the system is stable according to the 
Raylcigh and Schwarzchild criteria. 

An important feature in this context is that the oscillation 
spectrum is dense. This is because the local dispersion 
relation depends only on the wavenumber ratio Kr/K^. 
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As in a local analysis both \Kr and XK^ oo these may 
be adjusted to make the ratio take on any finite value and 
thus produce a continuous or everywhere dense spectrum 
(see eg. Papaloizou & Pringlc 1982; Lin ct al 1993 for a 
discussion of this point in a global context). 
This property of the spectrum means that commensura- 
bility conditions such as a = with fl being the angu- 
lar frequency associated with one orbit around a stream- 
line, required for parametric instability can always be ad- 
equately satisfied. 

For example (64) gives for incrtial modes {v^ = 0) and 
Keplerian rotation, tj^ = Vt^Kl/K"^. Thus when K^/Kr = 
l/y/3, a = ±n/2. 

Suppose, the radial and vertical wavenumber components 
have to be integral multiples of some smallest values 

{Kro, Kzq) respectively. Then because wc can let both 
wavenumber components become arbitrarily large, and be- 
cause of the ability to approximate irrationals arbitrarily 
closely with rationals, we can make \Kz/Kr — \/3| become 
as small as we like. 



Stobility Diogram 



Fig. 1. The stability diagram for equations (51) 

plied to the model with streamlines which are 



lipses in the (-jj^ 



and (52) ap- 
Keplerian el- 

e) plane. Instability occurs in the region of 



the plane between the two curves plotted. 



6. A steady state model with streamlines being 
Keplerian ellipses 

We now consider a more interesting non axisymmetric 
model for which the streamlines correspond to Keplerian 
ellipses with the same eccentricity circulating about a cen- 
tral point mass. This is described in section 3. In this 
particular model the streamlines form a pattern that is 
stationary in the inertial frame. However, in more gen- 
eral cases the pattern may precess slowly and so appear 
stationary only in a slowly rotating frame. But this ef- 
fect should make little difference. We comment also that 
because we are considering local stability, where pertur- 
bations arc localized along streamlines, these streamlines 
are only required to fill a small local volume of the fluid 
independently of exterior conditions. 
In this model fluid state variables vary as a fluid clement 
circulates around a streamline. Thus periodic forcing of a 
perturbed fluid element occurs and parametric instability 
becomes possible. 

The velocity field for the streamlines was specified through 
equations (7) - (8) in section 3. To proceed we need to 
specify two invariants for the streamlines. Because for the 
cylindrical disk model adopted, the motion is in horizontal 
planes so one of these may be taken to be z. The other 
may be taken to be a(l — e^) = r(l -|- cos(p)) which is 
proportional to the square of the angular momentum. 
Using these we adopt for the phase function 
S = Krr{l + e cos if) + K^z with Kr and again being 
constants. 

The wavenumber is then given by 

k = WS = {XKr{l + ecosf^), —XKrCsin^p, XK^). 

In this case 



r(l -|- ecos If) 



Moximum Growth Rote 



Fig. 2. Maximum growth rate in units of -y/ GM, / plotted 
as a function of eccentricity for solutions of equations (51) and 
(52) applied to the model with streamlines which are Keplerian 
ellipses. 



so that equation (37) gives 

e sin (fr+(l + e cos (f)(p 



bi = 



-y/(l + -f 2ecos(p) 



(66) 



Prom equation (38) we obtain 

Kz{l + ecos<^)f — Kzesmip<p 



b2 = 



N-s/{l + e'^ + 2ecosip) 



esmip,\ ,0 ,(65) 



Kr^y{l + e'^ + 2e cos ip)ii 



N 



(67) 
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where k is the unit vector in the vertical direction and 7. Solution of the governing equations 

N = J K?:(l + + 2ecos(p) + Kl is the normalizing fac- ^ ^ ^ ^ • ^- tr^ cr,^ u 

^ . J . . / • i We nrst note that the governine; equations 51-52 can be 

tor required to ensure that 02 IS a unit vector. -ij. • i. i. t ■i.-u ■ ■^ i i. 



written in a compact vector form with similar structure 
to the Lynden-Bell & Ostriker (1967) formulation of the 
general stability problem (see equations (18) and (24)): 



Given these unit vectors the forms of equations (51) and 
(52) can again be found by direct calculation. As the de- 
tails are somewhat tedious they are relegated to appendix 
1. y-2(7oe3 xy = E-y. (77) 

The coefficients (Ei^A = 1,2,7 = 1,2, Cn) contained in j.i ^ j- • i ^ t / o^ i i-i 

, ^. '-' . /.^N , r , r, Here the two dimensional vector y = (a, p) and the vec- 

the linear governing equations (51) and (52) tor a and p , . / a • \t ^ j- • i j. • 

° ^;or ea X y = (— p,a) . Ihe two dimensional matrix E 

are given y elements En, E12, E21, E22 as defined in section 4.4 

GM^ (l + 2ecosy + e^(l-3sinV)) above. 
■^11 — 1 ^ g2 _|_ 2e cos tp When a system has stream lines perturbed slightly from 

circles and the unperturbed circular streamline case has 
GAf* (l + ecos(/3)^ an underlying spectrum of oscillation modes with a dense 

^3 (1 + + 2ecos(p)^ ' spectrum, it is possible to show analytically that the sys- 

tem has a generic parametric instability. This is done in 
^ _ GM^ 3Kzesm(p{l + ecostfi) ^ ^ DCq appendix 2. 

N (1 + e'^ + 2e cos if) 12 Dt However, we have found that the most direct method of 

assessing the stability of a general system, with stream- 
^ _ GM^ S-ftT^e sin (/^(l + ecosLp) ^ ^ b + (70) li^^^s differing significantly from circles, governed by equa- 
N [1 + + 2e cos ip) ^ ^ Dt tion (77) is to integrate the system as an initial value 

problem starting from general initial conditions and to ex- 
aiid amine whether the solutions manifest exponential growth 

GM* K^{2 + 4ecosyj + 6^(2 - Ssin^ (p)) ■■ (se below). In this way, which is effectively exploiting the 

E22 = - _(_ e2 -)- 2ecos(^) b2b2.(71) Lagrangian nature of the system and solving for condi- 

tions in the neighbourhood of a fluid element as it circu- 
In addition lates around streamline, problems with boundary condi- 

I qj^ q \3/2 tions are avoided. 

Go = bi • 1)2 = \ — '"Q^ ^j2) However, before presenting an example we briefly com- 

V r N{1 + e +2e cos if) ^^^^ ^^^^ ^^le system (77) can be dealt with in an equiva- 

a,nd Isnt Eulerian formulation in which an eigenvalue problem 

of the same type as that in the Lynden-Bell & Ostriker 
-P^o = GM^ K^{l + e cos y)e sin formalism is recovered. To obtain the Eulerian formula- 

Dt N{1 + + 2ecosy) tion the convective derivatives in (77) are replaced as in 

/1 , ■:> , ■> \\ (15)- Thus 

(1 + 3c- + Ac VPS p) R;.il + ci-o>,p) \ ^ ' 

(l-he2-h2ecosv') J ^ > Dy ^ dy 



with 



As our system has spatial variation only around stream- 



DN _ IGM* K^esm(p^/l + ecos(p lines, we can write v Vy = \\\dy/ds, with s measuring 

~ ~y J.3 77 j^2 ^ j^2n _[_ g2 ^ 2ecosy) " length along a streamline. The problem is thus a 

one dimensional one with periodic boundary conditions in 
we also have s. One now adopts a time dependence through a factor 

expiiat) and obtains an eigenvalue problem for a of the 
^ ^ _ GM.(l + ecosy>) ^ fo^^ (24) namely 



{1 + e cos prK^ + K^e^ sin^ (p{l + e^ + 2e cos ip) 



-a'y + aB{y)+C{y)^0, (79) 
where each of the operators B, C is self adjoint with the 



jY2Q _|_ g2 I 2ecos<z')2 ' 

^ inner product for general vectors yi,y2 defined as 

Kfe^ sin^ p 



N'^{Kl + A'2(l e2 2ecos(p)) 



(75) (y2,yi) = j py*2-yiAd8, (80) 



together with where A is the stream tube area function (see section 3.2) 

and the integral is taken around the streamline of interest. 
/ GM* DN Kz{l + ecosp)^^^ This formulation allows us to note some of the results of 

b2 • bi = Y Jjt N'^{1 + e'^ + 2e cos p) ' Lynden-Bell and Ostriker (1967) such as that eigenvalues 

occur in quartets with complex conjugates and with op- 

which enables all the coefficients to be calculated. posite signs. 
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7.1. Numerical solution for Kepler Ellipses 

We have solved the governing equations (see equation (77) 
in section 7 numerically as an initial value problem. To 
do this we specify the ratio K^/Kr and integrate the 
system forward in time from arbitrary initial conditions. 
Instability then shows up as exponential growth, the rate 
may be easily estimated in practice. This procedure is very 
much simpler than attempting to solve a matrix eigenvalue 
problem. 

In the limit of small eccentricity the problem can be solved 
analytically using a standard approach to parametric in- 
stability problems. This is done in appendix 2. In this case 
the maximum instability occurs when K^jK^ = l/\/3. At 
this wavenumber ratio, the maximum growth rate is 

A = 3ef2/16. (81) 

The calculation is a Lagrangian one so viewed from the 

point of view of a fluid element circulating around a 
streamline, the instability occurs through the interac- 
tion of two inertial modes with frequencies ±Cl/2. In an 
Eulerian sense, if the pattern is stationary in the inertial 
frame, the frequencies as viewed there are —mifl + Q/2 
and —1712^1 — where (mi, 7712) are the associated az- 
imuthal mode numbers. These frequencies must be equal 
in the stationary frame because there is only one eigenfre- 
quency a as viewed from there as considered above. This 
means that mi — m2 = 1. Otherwise the rrii are not con- 
strained. 

Numerical solution of the governing equations gives the 

stability diagram plotted in figure 1 . This shows the region 
of the {Kz/Kr, e) plane with instability for < e < 0.7. 
When e — > 0, the numerical and analytic results agree. 
The region of instability collapses onto K^JK,. = l/\/3 
and the maximum growth rate agrees with equation (81). 
The band of instability increases with e. 
The Maximum growth rate in units of ^ GM^, / is plot- 
ted as a function of eccentricity in figure 2. This indicates 
that the maximum growth rate for a given e, is such that, 
to a good approximation, as in the small e limit, A oc e for 
e < 0.7. 

8. Discussion 

In this paper we have formulated a general local stabil- 
ity analysis applicable to astrophysical flows that appear 
steady in an appropriate uniformly rotating coordinate 
system. Such a flow might correspond to a stellar mode 
of oscillation that is a travelling wave in the azimuthal di- 
rection or an accretion disk with a free eccentricity. In the 
simple example we considered the disk streamlines were 
Keplerian ellipses. 

The general stability analysis considered perturbations lo- 
calized in the neighbourhood of unperturbed fluid stream- 
lines. The spatial localization was assumed to tend to zero 
as the wavenumber tended to infinity but more slowly than 
the inverse magnitude of the wavenumber itself allowing 
the adoption of wave packets. It is possible to consider 



the limit of infinite wave number taken at finite oscilla- 
tion eigenfrequency and obtain spectral properties from 
the local theory that apply to the global problem. The 
physical reason for this is that the degrees of freedom un- 
der consideration correspond to inertial or gravity modes 
whose frequency depends only on the ratio of wavenum- 
ber components. The group velocity thus approaches zero 
with localization in a way that ensures that localization 
may be maintained for arbitrarily long times. This justifies 
the local procedure. 

We found that the local stability of fluid on a particular 
streamline was governed by a pair of second order differen- 
tial equations, which in a Lagrangian description, evolved 
the perturbation associated with a particular fluid element 
as a function of time. This reduces the stability problem to 
an initial value problem rather than an equivalent eigen- 
value problem obtained in an Eulerian representation. 
When conditions do not vary around streamlines one ob- 
tains standard oscillation modes. However, in the more 
general case when conditions vary around streamlines the 
ordinary differential equations are linear with periodic co- 
efficients and the system was shown to become gencrically 
parametrically unstable. The dense spectra of inertial or 
gravity modes ensures that resonance conditions may al- 
ways be satisfied. 

We applied our considerations to study the stability of 
fluid for which the streamlines were in Keplerian ellipses 
with eccentricity up to 0.7. We showed that for small e, 
the maximum instability occurs when K^/Kr = 1/ a/S and 
involves, in the Eulerian description, two modes with az- 
imuthal mode number differing by unity in magnitude. 
Then the growth rate is 

A = 3efi/16. (82) 

For large eccentricities, to a good approximation, this 
gives the maximum growth rate as a function of eccen- 
tricity. However, the band of unstable wavenumber ratios 
shifts (see figure 1). 

Finally the considerations in this paper are based on linear 
theory and so they cannot determine the outcome of the 
instability. Results based on three dimensional numerical 
simulations for accretion disks with elliptical streamlines 
will be presented in a following publication. 
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APPENDIX 1 



Keplerian streamline model 

Wc here consider the derivation of the coefficients in the 
linear governing equations for a and /? for the model with 
streamlines corresponding to Keplerian ellipses described 
in section (6). The velocity field is given by equations (7) 
- (8) and the two invariants on streamlines are taken to 
be z and o(l — e^)=r(l + cos<^)) being proportional to 
the square of the angular momentum. 
We here note that the well known orbit equation given 
above may be used to eliminate r in the expressions given 
below so that the coefficients only depend on 99 and quan- 
tities constant on a streamline. Given the dependence of if 
on time, they may also be expressed as functions of time 
and quantities constant on a streamline. 
The phase function is 

S = Krr{l + ecostp) + K^z with Kr and being con- 
stants. 

The wavenumber is given by 

k = AV5 = {XKr{l + ecos(p), —XKresinip, XK^). 

In this case 



Prom equation (38) we obtain 



I y r(l + ecosif) 
so that equation (37) gives 



Kz{l + e cos — A'_f:: aunpip 
7V-v/(l + e2 + 2ecosv?) 



Kr^J{l + C- + 2( (■OS^)k 

N ' 



(.85) 



where k is the unit vector in the vertical direction and 
N = ^ (l + + 2e cos <^) + K1 is the normalizing fac- 
tor required to ensure that \>i is a unit vector. 
Given these unit vectors the forms of ccjuations (51) and 
(52) can again be found by direct calculation. 
For this purpose it is useful to note for this model that 
the operator 



D 
Dt 



vV 



d 



a/ (1 + ecosip)- — h 



re sm ip 



(.86) 



and to recall that 

dr/df = (p and dCpjd^ = — f, while the unit vectors have 

zero derivative with respect to r. 

Applying this operator to equation (37) we obtain 



GAf, 



(1 + ( cos 



~Dr~V"r3 (l-|-e2-|-2ecos<^)3/2 
(—(1 + ecos<^)f + esinf^yj) 
and applying it to equation (38) we get 

D{N\i^) IGM^ / isT^ (1 + e cos ^)3/2e sin 



(.87) 



Dt 



(l-|-e2 + 2ecos(^)3/2 



^(l + e2 + 2ecos(^)3/2'^ 

^ V(l + e2 + 2ecos(p) J ' 
with 

DN IGM* isT^esinvsyrTe^os^ 



Dt V r3 + (1 + e2 + 2e cos cp) ' 

Prom the above relations we obtain 



GM*(1 + ecosw) „\ , 
esm^,V ^,0 ,(.83) 



Go = bi • b2 = 



GM* K^{1 + e cos ipf/^ 
N{l + e'^ + 2ecosip) 

_ GM^ (l + ecosy;)3 

(1 + e2 + 2ecos<^)2 ' 



(.88) 

(.89) 

(.90) 
(.91) 



bi = 



e sm u:r 



cos ip)(p 



y^(rTe2"+2ecos^ 



(.84) 



and 
b2b2 



• 2 GM«(1 + ecosv?) 
bo = : 
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(1 + e cos ^pYKi + K;e^ sin"^ (^(1 + + 2e cos ip) 
iV2(i + e2 + 2ecos</?)2 

sin2 LP 



N^{K^ + K'^{l + e^ + 2ecos^)) 



in the linear governing equations (51) and (52) for a and 
p. These can be constructed using equations (46-54) of 
section 4. We obtain: 

GM* (1 + 2e cos + 6^(1 - 3 sin^ (^)) 

-'^ii ~ ^ T"; — 

(92) r"* 1 + + 2e cos ip 



GM* [l + ecosipf 

We now use the general results for any pair of unit vectors — ^T\~, — i , o \2' V^'^^) 

,^ , J' (1 + + 2ecos(p)^ 

obtained from (53-54) that ^ 

GAf* SKrC sin (^(1 + e cost/?) • • DCq 
b2-bi=-bi.b2 + -^ (.93) ^12= ^3 ^(i+e2 + 2ecosv) +bi •b2 - -^^,(-102) 

and GM, 3K,e sin (^(1 + e cosy') ^/ K^^<^o^1nQ^ 

K K K K ,,,, ^^^-^ iV(l+e2 + 2ecos^) +^1 + ^(-103) 

bi • b2 = -b2 • bi — (.94) 



Dt 



and 



to obtain b2 • bi and bi • b2 which occur in the equations r^c/r. a ores ^ ■ i w 

A a GM* isT^ (2 + 4e cos + 6^(2- 3 smV)) , i: , 

governing a and /?. E22 = —3 — — ,210 ^ ^-b2b2.(.104) 

We first apply the convective derivative operator to (.90) iV (i + e + .^e cos (pj 

to obtain 

DCq GMt, Kz(l + e cos ip)e sin 
Dt ~ N^l + e^ + 2ecosLp) ^ 

(l + 3e2-|-4ecos(^) K'^{1 + ecosip)\ 

(.95) 



(l-|-e2 + 2ecos(/9) N'^ 
and use (.87 - .89) to obtain 



h h = fGM:^^ j^,(l + ecosy)3/2 
' ^ ^ V r3 £)t 7V2(i + e2 + 2ecos<^) ^ 

which enables these quantities to bo found directly. 

For this model, because the velocity field is Keplerian, we 

have the very useful result that 

^•^^ = -^r. (.97) 

Using this together with ^ given by equation (39) we ob- 
tain 

u V7/ V7 ^GM, 3Js:^e sin (^(1 + e cos (p) 

• • • = jV(l + e2 + 2ecosy) 

GM* 1 + 2e cos + e2(l - 3 sin2 (^) 

—a — 5 ^5 — , (.98) 

r'' 1 + + 2e cos y 

and 

GM* ZK^e sin (^(1 -|- e cos <^) 



b2 • (^ • V(v • Vv)) = a- 



N{l + e'^ + 2ecosip) 



CM, Kl{2 + 4e cos ^ + 6^(2 - 3 sin^ ^)) 
^ r3 iV2(i + g2 + 2ecos(^) ^" ^ 

In addition for this particular model with uniform pres- 
sure, we automatically have for i = 1, 2, j = 1, 2, 

We have now assembled all the quantities needed to calcu- 
late the coefficients {Co,Eij,i = 1,2, j = 1,2) contained 
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APPENDIX 2 

Parametric instability for nearly circular stream- 
lines 

We here consider the generic parametric instabihty that 
occurs for streamlines near to circles for systems with iner- 
tia/gravity modes. But note that as indicated for the spe- 
cific example of Keplerian ellipses instability may persist 
for significant departures of the streamlines from circles. 
We start with the governing equations in vector form (see 
equation (77) in section 7). These are written as: 



y - 2Coe3 x y = E • y, 



(.105) 



where the two dimensional vector y-^ = {a,f3) and the 
vector es X y = {—$, oi)^ . The two dimensional matrix E 
has elements Eii,Ei2,E2i,E22 as defined in section 4.4 
above. 

For circular streamlines E = (E)^ and Cq = Co,c are 
constant, while for small departures from circularity we 
can expand to first order in some small parameter e. 
When the streamlines are Keplerian ellipses e = e. We 
thus write Cq = Co.c + £Cq,i and E = Ec + e£. Here 
Co,i and £ are periodic functions of time where the pe- 
riod is just the time it takes a fluid element to orbit 
around a streamline. Without loss of generality we may 
assume that the time average of these coefficients over 
one period is zero. For the Keplerian ellipse streamlines, 
from the results given in appendix 1, Co,c = ^K^/K, 
and the only non zero element of (E)^ is (£^22)0 = 
3{flKz/K)'^. Here Q, = ^JGM^|a^ where a is the semi 
major axis of the streamline of interest. Similarly, mak- 
ing the equivalence e = e, Co,i = il^K^j K^" co&ip while 
the elements of £, are given by £\\ = — O^cost^, £12 = 
n2(i;f^/ii'3)sin<^(4if2 + 2if2)^ = 2n'^{K,/K)smip 

and £22 = 2n^{K^/K^) cos<f{AK^ + Kj). 
Note that we may write to this order = Q.t making 
the coefficients that are proportional to the eccentricity 
vary harmonically with time. The period is simply the or- 
bital period on the streamline. We now write the governing 
equations in the form 

y - 2Co,ce3 x y - (E)^ • y = 2eCo,ie3 x y + • y (.106) 

and treat the right hand side using perturbation theory. 
We first remark that when e = 0, there is a simple normal 
mode solution appropriate to the time independent case 
which may be written in the form 



y = yoexp(ic7t), 



where yo is a constant vector and cr is a normal mode 
frequency which satisfies 



-CT^yo - 2if7Co,ce3 x yo - (E)^ • yo = 0. 



(.108) 



In line with the example presented in section 5, we expect 
this mode to correspond to an inertial/gravity mode. The 
spectrum is dense and so we may always expect to be able 



to choose (T to be close to 0/2 which is the lowest order 
commensurability condition for parametric resonance (for 
the example of section 5 the commensurability is exact for 

We now move on to consider the case of e non zero but 

small and set a = 0/2 + 5a, where Sa = 0(eO). 

We now seek a solution of equation (.106) in the form 



y = A{t)yQ exp(iOt/2) + cc. + Y^ e"y„ 



(.109) 



Tl=l 



Here A{t) is a slowly varying amplitude function to be 
determined and here and below cc denotes the addition 
of the complex conjugate of the first expression ensuring 
that as required y is real. The terms e"y„, assimied real 
are additions that are of first and higher order in e. 
By stating that A is slowly varying we mean that \A/A\ = 
0(fO) and \A/A\ = 0(eO)2. 

Given this we substitute equation (.109) into equation 
(.106) and drop expressions that are clearly second or 
higher order in e to get an equation for A. This reads 

(iOyo - 2Co,ce3 x yo) {^A — i5<jA^ exp(iOi/2) + cc = 

e (iOCo^es x yo + £yo) A exp(iOt/2) -|- cc 

-e (yi - 2Co,ce3 x yi - (E)^ • yi) (.110) 

We now multiply by )exp(— zOt/2) and take a time aver- 
age over many orbital periods, assuming A varies slowly 
enough that it can be factored out, to obtain 

(iOyo - 2Co,ce3 x yo) {A - iScjA^ = 

(e (-iOCo,ie3 x y^ + £yl) exp(-iOt)) A* 

-e((yi - 2Co,ce3 x yi - (E)^ • yi) exp(-iOt/2))(.lll) 

where the angular brackets denote the time average over 
many orbital periods. Because Co,i and £ vary harmoni- 
cally with angular frequency O, it is anticipated that the 
time averages retained above are not zero. 
Setting A = ^oexp(At) in equation (.111), A being the 
growth rate, taking the scalar product with yQ, which as 
can be shown by doing the time integration by parts re- 
duces the final term involving yi to be O(e^) and finally 
taking the modulus of the result gives an expression for 
correct to O(e^) in the form 



(.107) A^ = 



Keyg-gyg exp(-»Ot))|' 
|iO|yo|2 - 2Co,cyo ' ea x yo| 



(.112) 



Thus parametric instability occurs (A is real) at a maxi- 
mum rate when the commensurability is exact or a = 0/2 
and for fixed e the frequency half-width associated with 
instability corresponds to this maximum growth rate in 
magnitude. The growth rate itself is oc e. 
It is easy to evaluate equation (.112) for the Keplerian 
streamline model of section 6. We begin by noting that 
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yo can be found from the solution for the local oscillation In this limit ^ is of the form (39) and is obtained from 
mode given for the circular streamline case in section 5 the eigenvalue problem (79), as is a, which then belongs 
(noting that O' = — 3il/(2r) for the Keplerian case). Prom to the spectrum of C 

that yI = {2iK^n,Ka). 

It is then straightforward to evaluate A given y, Co,c and 
£ as specified above. One obtains 



giving a peak growth rate as a function of eccentricity 
given by Xmax = SefJ/lG. The instability band width is 
then given by 6(j = Xmax- 

APPENDIX 3 

Relation to the Spectrum of the global problem 

We have studied the stability of fluid displacements that 
are localized around the vicinity of fluid streamlines and 
circulate around with the fluid. When conditions vary pe- 
riodically around streamlines, parametric excitation be- 
comes a possibility aided by an underlying mode spec- 
trum that is dense so that resonance conditions may al- 
ways be satisfied. The analysis was done in a local limit 
of high wavenumber in directions orthogonal to stream- 
lines. However, it is possible to relate the local spectrum 
so found to that of the global problem. A physical reason 
for this is that it is possible to take a local limit where 
wavcnumbers tend to infinity but mode oscillation fre- 
quencies remain finite. This naturally results in a limit 
in which the group velocity vanishes and hence localiza- 
tion becomes possible. We here illustrate the connection 
between the local and global problems. 
The equation satisfied by the Lagrangian displacement ^ is 
equation(18) which after use of equation (19) and replace- 
ment of the convective derivatives through equations(22) 
and (23) yields an operator equation of the form 



where we have included the <j dependence to emphasize 

that we have an eigenvalue problem for a and that the de- 
termined eigenvalues constitute the spectrum of the sys- 
tem. We adopt the inner product defined through (25) to 
define the norm |^| = |(^,^)|. 

Then following Papaloizou & Pringle (1982) and Terquem 
& Papaloizou (1996), we remark that the spectrum coin- 
cides with the values of a for which the operator C has 
an unbounded inverse. Further if one can find any ^ for 
which |£(^, C7)|^/|^|^ can be made as small as we please 
then C will have an unbounded inverse and so a belongs 
to the spectrum. 

In fact displacements of the form (26) which are localized 

on streamlines over a length scale that tends to zero as 
a fractional power of the wave number parameter A~^, 
which itself scales as the inverse wavenumber as that ap- 
proaches zero, can be used to construct a set of normalized 
^ which lead to |£(^, cr)p as (see Papaloizou 

& Pringle 1982; Terquem & Papaloizou 1996 for other ex- 
amples). 



A^ = 96^07256 -((5a2), 



(.113) 



£(C,a)=0, 



(.114) 



